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ABSTRACT 

Using the Yunnan evolutionary population synthesis (EPS) models with and without 
binary interactions, we present the luminosity of Ha recombination line (Lua), the 
luminosity of [OIIJA3727 forbidden-line doublet (L[oii]), the ultraviolet (UV) fluxes 
at 1500 and 2800 A (Li,uv) and far-infrared flux (Lfir) for Burst. SO, Sa-Sd and Irr 
galaxies, and present the calibrations of star formation rate (SFR) in terms of these 
diagnostics. 

By comparison, we find that binary interactions lower the SFR.vs.Lhq and 
SFR.vs.L[oii] conversion factors by ~0.2dex. The main reason is that binary interac- 
tions raise the UV flux (shortward of the Lyman limit) of the stellar population (SP) 
in the age range 6.7< logt/yr <8.4 and thus more ionizing photons are present in the 
nebula. Moreover, binary interactions do not significantly vary the calibrations of SFR 
in terms of Li.uv- This is because binary interactions raise the flux at 1500 A of the 
SP in the range 8.75 < logi/yr < 9.2 and the maximal difference is about 1 dex. In 
addition, binary interactions have little effect on the flux at 2800 A. At last, the cali- 
bration of SFR from Lfir is almost unaffected by binary interactions. This is caused 
by the fact that binary interactions almost do not affect the bolometric magnitudes 
of SPs. 

We also discuss the effects of initial mass function (IMF), gas- recycle assumption 
and EPS models (including gisselqs, boos, starburstqq, PopSTAR and pegase mod- 
els) on these SFR calibrations. Comparing the results by using Salpeter (S55) IMF 
with those by using Miller & Scalo (MS79) IMF, we find that the SFR.vs.Lna and 
SFR.vs.L[oii] conversion factors by using S55 IMF are greater by 0.4 and 0.2 dex than 
those by using MS79 IMF for the Yunnan models with and without binary interactions, 
respectively. The SFR.vs.Li^uv and SFR. vs. Lfir conversion factors by using S55 IMF 
are larger by an amount of 0.2 dex than the corresponding ones by using MS79 IMF. 
The inclusion of gas-recycle assumption only lowers these SFR calibrations at faint 
SFR. Moreover, comparing the results when using different EPS models, we find that 
the differences in the SFR.vs.Lhq and SFR.vs.L[oii] conversion factors reach 0.7 and 
0.9 dex, the difference in the SFR. vs. Lfir conversion factor reaches 0.4 and 0.8 dex, 
and the differences in the SFR.vs.Lj uv conversion factors reach 0.3 and 0.2 dex when 
using S55 and NON-S55 IMF (including ChaOS, KOI, K93' and MS79 IMFs, partly 
caused by the difference in the IMF), respectively. At last, we give the conversion 
coefficients between SFR and these diagnostics for all models. 

Key words: binaries: general - galaxies: fundamental parameters - galaxies: general 



1 INTRODUCTION 

One of the most recognizable features of galaxies along the 
Hubble sequence (loose definition, including not only mor- 
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phological type but also gas content, mass, bar structure 
and dynamical environment) is the wide range in young 
stellar content and star formation activity. Understanding 
its physical nature and origin of the variation in stellar con- 
tent are fundam ental to understand evolution of galaxies 
(|Kennicuttlll998l . hereafter K98). Star formation rate (SFR) 
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can be used to compare with those distant galaxies at cosmo- 
logical lookback times, and extrapolate the future timescales 
for star formation in galaxies by combining: with HI ari d 
CO measurements (jKennicutt. Tamblvn fc Congdonlll994 ). 
Moreover, star formation activity is usually correlated with 
cold gas and stars in galaxies: stars continuously produce 
mass, energy and metals during their evolution processes, 
and return them to galactic medium (gas), affecting the sta- 
tus of the next generation of stars. SFR carries the informa- 
tion on the evolution of galaxies. Therefore, it is important 
to determine SFR and its variation with Hubble type (loose 
definition) and environment, which can help us to under- 
stand the evolution of galaxies. 

The commonly used SFR tracers include the flux of Ha 
nebular recombination line, the ultraviolet (UV) continuum 
flux, the flux of [OH] A 3727 forbidden- line doublet and far 
infra-red (FIR) continuum flux (K98). These SFR tracers 
are more or less correlated with the UV passband. 



• First, the UV flux is directly tied to the photospheric 
emission of the young stellar population (SP). 

• Second, the integrated luminosity of galaxy shortward 
of the Lyman limit (Far-UV) can ionize the hydrogen in the 
nebula and produce the recombination lines (such as Hq, 
H/?, and so on). Thus the luminosities of these lines can be 
used to trace SFR. 

• Third, the luminosity of the strong [OII]A3727 
forbidden-line doublet is often empirically obtained through 
the Ha luminosity, although it is not coupled to the ioniz- 
ing luminosity and the excitation of this line is sensitive to 
abundance and the ionization state of the gas. 

• Finally, the last SFR trace, the FIR luminosity, is also 
correlated with the UV passband. The interstellar dust can 
absorb the bolometric luminosity of galaxy and re-emit it 
in the thermal IR passband. The absorption cross section of 
the dust peaks in the UV passband, and, since the UV flux 
is considered as a tracer of young SP, the FIR luminosity 
can also diagnose SFR. 



Furthermore, the last three diagnostics and the corre- 
sponding calibrations of SFR are correlated with the UV 
flux. Since in our previous studies, we found that the inclu- 
sion of binary interactions in evolutionary population syn- 
thesis (EPS) models can raise the UV flux by ^ 2-3 magni- 
tudes for SP at an age of ~ 1 Gyr jZhang et al.ll2004 , ,200^ , 
in this study we will discuss the effect of binary interactions 
on these calibrations. 

The outline of the paper is as follows. In Section 2 we 
describe the used EPS models and algorithm. In Section 3 
we overview the previous results about SFR calibrations, 
and the advantages and disadvantages of these SFR tracers. 
In Section 4 we give the effect of binary interactions on these 
SFR calibrations. In Section 5 we discuss the effects of initial 
mass function (IMF), gas-recycle assumption and the EPS 
models on these SFR calibrations, and give the conversion 
coefflcients between SFR and these tracers for all models. 
Finally we present a summary and conclusions in Section 6. 



2 MODELS AND ALGORITHM 
2.1 Spectral synthesis models 

First, we use the Yunnan EPS models, w hich have been bui lt 
by Zhang and her colleagues since 2002 (jZhang et al.ll2002l ). 
The main characteristic of Yunnan EPS models is the in- 
clusion of various binary interactions. Yunnan EPS models 
have given the results of SPs with and without binary inter- 
actions at seven metallicities (from 0.0001 to 0.03) and 90 
ages (from O.lMyr to 15Gyr). 

The Yunnan EPS models were built on the basis 
of the Cambridge stellar evolutionary tracks jEggletonl 
1 19711. 1 19721. 1 19731). BaSeL-2.0 ste l lar a tmosphere models 
( Leieune. Cuisinier fc Buseij Il997l . 1 19981 . hereafter LCB97, 
LCB98) and various initial distributions of stars. The Cam- 
bridge stellar evolutionary trac ks are obtained by using 
the rapid stellar evolution code (|Hurlev. Pols fc Toutll2000l : 
iHurlev. Tout fc Polsll2002l). which is based on the stellar evo- 
lutionarv tracks bv lPols et al.l (|l998h . 

In this work we use a set of standard Yunnan EPS models 
at solar metallicity. The description of standard models is 
as follows: (i) the initial mass of the primary Mi is given by 

^ ^ (1 - X)0'75 + 0.032(1 - X)0'25 ' ^ ' 

where X is a random variable uniformly distributed in 
the range [0, 1]. This exp r ession is the approximation by 
lEggleton. Fitchett fc 'IbutI ll 19891. hereaft^er EFT ) to the 
IMF {(f,{M) = diV/dM) of iMiUer fc Scald (|l979l . hereafter 
MS79): 

r M"^ '', 0.10 <M< 1.00 
(t>{M)^s^g oc I M-2■^ 1.00 < M < 10.0 (2) 
i M-^■^ 10.0 < M < 100. 

in which M is the stellar mass in units of Mq; (ii) the ini- 
tial masses of the component stars in a binary system are 
assumed to be correlated, and the initial mass of the sec- 
ondary is obtained from a uniform mass-ratio q distribution 
(|Mazeh et al.lll992l : IColdberg fc Mazeh|[l993 ) ; (iii) the dis- 
tribution of separation between two component stars takes 
the following form: 

in which a scp ~ 0.070, ap = lOR^.O! = 5.75 x lO^'Ra and 
m ^ 1.2 (|Han. Podsiadlowski fc EggletonI llQgsl ): and (iv) 
the eccentricity values follow a uniform distribution. 

In the standard Yunnan EPS models, approximately 50% 
of stellar systems are binary systems with orbital periods less 
than lOOyr. This fraction is a typical value for the Galaxy, 
resulting in ~ 10.1% of the binaries experiencing Roche lobe 
overflow during the past 13 Gyr (see Han et al. 1995). 

To investigate the effect of IMF on the results, we use 
another set of solar-metallicity Yunnan EPS models, which 
differs from the standard models by the initial mass distri- 
butions of the primary and secondary stars. In this set of 
models, (i) the initial mass of the primary is given by 

Mi^OMj^f'"- (4) 

This exp r ession is the approximation by EFT to the IMF of 
ISalpeteij (|l955l . hereafter S55) with a = -2.35, i.e.. 
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Table 1. Description of the used Yunnan, BC03, GISSEL98, PopSTAR, STARBURST99-v6.02 and PEGASE-v2 EPS models. The superscript 
'C means that it is available to choose, 'A' means that it is added by this work. 



name 


evolutionary LIB., Z, age range (yr) 


spectral LIB. 


IMF 


Mi,Mu 


OUTPUT 










(Mo) 


ISED/Q/Lii„es 


Yunnan 


Cambridge ,solar, 0.1M-15G 


BaSeL-2.0 


S55/MS79*! 


~ 0.1-100 


Y/Y/Y 


BC03 


Padova'^ , solar, 1.0M-20G 


BaSeL-3.1'^ 


S55/Cha03 


0.10-100 


Y/N/N 


GISSEL98 


Geneva*^ ,solar, 1.0M-20G 


combination*"^ 


MS79 


0.10-100 


Y/N/N 


PopSTAR 


Padova ,solar, O.lM-6. G 


BaSeL** 


SSS'/KOl^ 


0.15-100 


N/Y/Y*5 


STARBURST99-V6.02 


Padova AGB'^ ,solar, 1.0M-15G 


BaSeL'^ 


S55/K93'^ 


0.10-100 


Y/Y/Y*e 


PEGASE-V2 


Padova, 1.0M-20G 


BaSeL-2.0+CM 


S55/K93'C 


0.10-100 


Y/Y/Y*7 



*1 In the Yunnan models the IMF is for the primary in binary a system. 

*2 GISSEL98 models use several stellar evolutionary libraries, Geneva is the main one. 

*3 GISSEL9S models use the combination of several stellar spectral libraries. 

*4 PopSTAR models use several stellar spectral libraries, BaSeL is the main one. 

*5 PopSTAR models provide the number of ionizing photons Q(H, Hel, Hell, OI), the luminosities of emission-lines i(Ha,H/3)- the 2nd 
of a series, the luminosities of the other 18 emission-lines are provided as a function of Q(H) for HII region. 
*6 STARBURST99 models provide Q(H, Hel, Hell) and i(Ha,H^,P/3,BG) ■ 
*7 PEGASE models provide Q(H) and the luminosities of 61 spectral lines. 



0(M)s55 = M- 



(5) 



and (ii) the masses of the component stars in a binary system 
are assumed to be uncorrelated, i.e., the secondary mass is 
chosen independently from the same IMF as the primary. 

In Table. [T] we give the main characteristics of the used 
EPS models. The first column gives the name of each EPS 
model, the 2nd-5th columns give the used stellar evolution- 
ary tracks (including metallicity Z and the age range), the 
used stellar atmosphere models, the used IMF and the lower 
and upper mass limits (Mi and Mu), and the last column 
describes whether the integrated spectral energy distribu- 
tions (ISEDs), the number of ionizing photons Q and the 
luminosities of spectral lines I/unes are provided by the cor- 
responding models. 



2.1.1 Other spectral synthesis models 

To check the effect of spectral synthesis models on 
these SFR calibrations, we also use the GISSEL98 
(Gala xy Isochrone Synthesis Spect ral Evolution Li- 
brary , iBruzual fc Chariot! Il993l), BC03 jBruzual fc Chariot 
200s ) , PopSTAR (jMolla. Garcia- Vargas fc Bressan 
200g). STARBURST 9 9 v 6. 02 dLeitherer et al.1 llQOj 
Vazquez fc Leithererl 20051: Leithcrcr ct al. ' I2OIO!) and 



PEGASE y2 (jFioc fc Rocca-Volme ranec 1997, 1999) models. 
All these models do not include binary interactions, and 
the characteristics of them are also summarized in Table [T] 
Next we describe in more detail the main ingredients of 
these additional EPS models. 

a. GISSEL98 and BC03 models 

GISSEL98 and BC03 models were built by Bruzual fc 
Chariot in 1993 and 2003, respectively, and provided the 
results of SPs at six metallicities (from 0.0001 to 0.05) and 
221 ages (from 1 Myr to 20 Gyr) in tabular form. 

GISSEL98 m odels were based on the Geneva stellar evo- 
lutionary tracks (jMaeder fc Mevnet|[l989l . Il99ll ) mainly, the 
MS79 IMF with the lower and upper mass limits Mi = 0.1 
and A/u — 100. Mq (see Eq. [2]), and the combination of sev- 



eral stellar spectral libraries (we refer the interested reader 
to their paper for details). 

For the BC03 models, they used two sets of stellar evo- 
lutionary tracks (padova-1994 & Padova-2000), two forms 
of IMF [S55 with a = - 2.35 fc Chabrier (200 3, here- 
afte r Cha03 ) ], hig h- [STELIB (|Le Borgne et al.ll2003l ) fc Pick- 
les jPicklesI 1 19981 )] and low-resolution [BaSeL-3.l, LCB97, 
LCB98] stellar spectral libraries. The Cha03 IMF is as fol- 
lows: 



fCiM- 
l C2M- 



-(logM-logMc)^ 



exp 
2.3 



(6) 



M < 1.0 
M > 1.0 

in which Mc = 0.08 Mq, a ^ 0.69 and M is the stellar mass 
in units of Mq. For both IMFs, the lower and upper mass 
limits are 0.1 and lOO.M©. 

For the GISSEL98 models, we use the set of ISEDs for 
solar-metallicity SPs. For the BC03 models we use the set 
of ISEDs of solar-metallicity SPs, which are derived by us- 
ing the Padova-1994 stellar evolutionary tracks and BaSeL-3.l 
stellar spectral library. 

b. PopSTAR models 

The PopSTAR model s w ere built by 
iMolla. Garcia- Vargas fc BressanI (|2009|). In their models, 
they used a revision of the Padova (|Bressan. Granato fc Silva 
Il998h isochrones used in iGarcia- Vargas. Molla fc Bressan 
()l998h . and the BaSeL stellar atmosphere models (LCB97, 
LCB98) mainly. All the results [the number of ionizing pho- 
tons Q(H, Hel, Hell, 01) and the emission-line luminosities 
L{Ha,Hp)] of SPs are provided in tabular form, including 
six metallicities, six IMFs and 93 ages (from O.lMyr to 
logt/yr =9.78). We choose the two solar-metallicity sets, 
corresponding to the following IMFs: 

• the S55 IMF with a = —2.35 and mass range of 0.15 — 
100 M0 (because the Mi is different from that of the other 
models, we caU S55' IMF in Table [U; 

• the Kroupa IMF (2001, hereafter KOI) 

r CiM-°-^°, 0.01 < M < 0.08 
0(A/)koi = \ C2M-i-3°, 0.08 < M < 0.50 . (7) 
I C3M-2-3°, 0.50 < M < 100. 
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in which M is the stellar mass in units of M0 . 

Moreover in t he 2nd paper of a series 
iMartm-Manion et alj |2010| ). the luminosities of other 
18 emission-lines are given. 

c. STARBURST99--VG.02 code 

STARBURSTQ^ is a web based software and data pack- 
age designed to model spectrophotometric and related prop- 
erties of star-forming galaxies. It was developed by re- 
searchers in Space Telescope Science Institute. We use the 
6.02 ver sion. The description o f the model input physics is 
give n by Leitherer et all (Il999l '). lvi,zauez fc Leitherer (120051 ) 
and iLeitherer et al. I (|2010D " 



• For the STARBVRST99-yQ.Q2 code, four sets of steUar 
evolutionary tracks are provided, each set including 5 metal- 
licities. We choose the set of solar- met allicity Padova AGB 
tracks (including thermally pulsing AGB stars, the 44th 
tracks) . 

• Five sets of stellar atmosphere libraries are provided, 
we choose the BaSeL stellar atmosphere models (LCB97, 
LCB98). 

• By default, the Kroupa IMF with two mass intervals 
(the exponent a = [—1.3, —2.3], the mass boundary Mcut = 
-I- [0.1 ,0.5, 100.] M0) is used. To be consistent with those 
of the other spectral synthesis models, we change this IMF 
form to (i) the S55 IMF with a = -2.35 and (ii) the Kroupa 
(1993, hereafter K93) IMF with three mass intervals: 

-1.3 



r CiM~'-\ 0.10 < M < 0.50 

';*(^)k93 = { ■^ 0.50 <M < 1.00 

I CsM-^-^ 1.00 <M < 100. 



(8) 



in which Ci = 0.035, C2 = 0.019, C3 = 0.019 and M is the 
stellar mass in units of Mq . Because all coefficients in Eq. [8] 
are set to 1 in this study (see also pegase), we call K93' 
IMF in Table [U 

• Two cases of star formation (const, without) are pro- 
vided to choose, we choose the case of fixed mass (i.e., with- 
out star formation). 

At last, we obtain the ISEDs, the number of ioniz- 
ing photons (3(H, Hel, Hell), the emission-line luminosities 
i(Hc«,H;8,p^,BG) of solar-metallicity SPs at 83 ages in the 
range lMyr-15Gyr at an interval of logf/yr = 0.05. 



evolutionary tracks (7 metallicities, Z=0.0001, 0.0004, 
0.004, 0.008, 0.02, 0.05 and 0.1), the combination of 
BaSeL-2.0 (LCB97, L CB98 for T^s < 50 000 K) with 
IClegg fc MiddlemassI l| 19871 . hereafter CM, for T^tt > 
50 000 K) stellar spectral libraries are used. 

• Among the nine IMFs provided, we choose the S55 IMF 
with Q = —2.35 and K93 IMF. The upper and lower mass 
hmits of both IMFs are set to 0.1 and 100 M© (by defauh, 
Mu = 12O.M0), respectively. Because the coefficients of K93 
IMF (see Eq. [8]) are set to 1 in the PEGASE code, we also 
call K93' IMF in Table [T] 

• Among the six forms of SFR provided, we choose the 
instantaneous burst, the const SFR and the exponentially 
decreasing form {SFR = p2 ■ exp(— r/pi)/pi, p2 = 1.0) with 
the timescales (pi) of 1, 2, 3, 5, 15 and 30 Gyr. They are used 
to built Burst, Irr, E-Sd types of galaxies, respectively. 

• For the other model input parameters, the default val- 
ues are used. Metallicity (mass fraction) of the interstellar 
medium (ISM) at t = is zero, the fraction of close binary 
systems is 0.5, and the mass fraction of sub-stellar objects 
formed is 0.0. And the default evolution processes are used, 
the evolution of stellar metallicity is consistent, stellar wind, 
infall, galactic winds and global extinction are neglected and 
nebular emission is considered. 

Using the above set of input parameters and physics, 
we obtain the number of ionizing photons Q(H) and the 
luminosity of recombination line Lhq for Burst, E, SO, Sa- 
Sd and Irr galaxies at 68 ages in the range of lMyr-20Gyr. 



2.1.2 Comments on the above mentioned models 

For PEGASE models we directly use their results because 
they consider the star formation and nebular emission. 
For the other models {GISSEL98, BC03, PopSTAR and STAR- 
BURST99), we need to generate the results of galaxies with 
different galaxy types by means of spectral synthesis mod- 
els. Because the ISEDs are not provided for PopSTAR mod- 
els, we could not give the ISEDs of galaxies with difi^erent 
galaxy types (therefore the UV and FIR continuum fluxes), 
we only could give the luminosities of the Ha recombina- 
tion line and the [OII]A3727 forbidden- line doublet, and the 
nebular emission continuum by using the number of ionizing 
photons Q(H) provided by them. 



d. PEGASE-V2 code 

PEGASE is a code to com pute the spectral evolution 
of galaxies, a nd was developed bv lFioc fc Rocca-Volmerangel 
l| 19971 . 1 19991 ). The evolution of the stars, gas and metals is 
followed for a law of star formation and a stellar IMF. The 
stellar evolutionary tracks extend from the main sequence to 
the white dwarf stage. The emission of the gas in HII regions 
is also taken into account. We use the 2.0 version. The main 
improvement in version 2 is the use of evolutionary tracks 
of different metallicities (from 10~* to 5Z0). The effect of 
extinction by dust is also modeled using a radiative transfer 
code. 

• For the PEGASE-v2 code, a set of Padova stellar 



^ http:/ /www. stsci.edu/scicncc/starburst99/ 
2 http://www2.iap.fr/users/fioc/PEGASE.html 



2.2 Construction of galaxies with different galaxy 

types 

The construction of galaxies with differ ent galaxy types ha s 
been described in our previous paper (|Zhang et al.1 bOld ). 
In brief, we use the BC03 software package to build them. 
Eight galaxy types (Burst, E, SO, Sa-Sd and Irr) are in- 
cluded, and they are built by a delta-form SFR, six expo- 
nentially decreasing SFRs with characteristic time-decays 
T = 1,2,3,5,10,15, and 30 Gyr and a constant-form SFR, 
respectively. The exponentially decreasing SFR is given by 

i>it) = [l + eMpG(t)]r-'exp(-t/r), (9) 

where r is the e-folding timescale, MpG(i) = [1 — 
exp(— i/r)] — Mstars — Mremnants is the mass of gas that has 
been processed into stars and then returned to the ISM at 
time t, Mstars and Mromnants are the masses of stars and 
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remnants at t, and e denotes the fraction of MpG(t) that 
can be recycled into new star formation. 



2.3 Computations of Lua, i[oii)? ^fir and nebular 
continuum 

In this part, we will describe the computations of the lu- 
minosity of the Ha recombination line Z/Hq , the luminosity 
of the [OII]A3727 forbidden-line double L[on], the FIR flux 
LpiR and the emission of nebular continuum -Fncb.A- 

Before the computation of Lhq, we must calculate the 
number of ionizing photons Q(H). During the computation 
of Q(H), one method assumes that it only comes from young 
(< lOMyr) and massive (> 20-25 M©) stars (such as K98, 
STARBURST99) ■ Another method obtains that number by 
integrating the photons below 912A (such as PopSTAR). We 
choose the second method, i.e.. 



Q(H) 



dv, 



(10) 



where Fi, is the stellar flux in Hz, 1/ is frequency and h is 
the Planck constant (=6.6262 x 10"^^ ergs). 

Then, we assume that all the star formation is traced 
by the ionized gas, although in the PEGASE code 70 per cent 
of the Lyman continuum photons computed from the spec- 
tral synthesis models are absorbed by the gas (i.e., ionize 
the gas) and the rest (30%) are absorbed by dust when con- 
sidering extinction. At last, we use Case B recombination 
at electron temperature re=10,000K and number density 
= lOOcm"^. The same set of the parameter values (Case 
B, Te and rie) is used by K98, STARUBURST99 and other 
studies. Under the above assumptions, the luminosity of Hq 
can be obtained by the following expression (popSTAR): 



(11) 



where ob is the recombination coefficient to the excited level 
in hydrogen, which d epends on the electronic temperature, 
j-B and Or are from [ Ferlandl ()l98d ). and the ratio a/P is 
taken from lOsterbrock 



The luminosity of [OIIJA3727, L[oiii, is often obtained 
by an empirical method. In this paper we obtain it by as- 
suming the ratio of the luminosity of [OII]A3727A to Ha, 
L[oii]/Lu_a = 0.45, as used by K98, although in the PEGASE 
code it equals to 3.01/2.915, and the discrepancy in the 
L[oii]/Lhq is large for the different types of galaxies. 

The FIR luminosity is mainly used to calibrate the 
SFR of starburst galaxies, for which the FIR luminosity 
is often assumed to equal to the bolometric luminosity, 
ipiR ~ isOL. 

At last, the emission of nebular continuum can be ob- 
tained by: 



p 

ne 



Q(H), 



(12) 



where c is the light velocity and F is the emission coeffi- 
cient for hydrogen and helium (He/H=0.1), which includes 
free-free and free-bound contributions and the emission co- 
efficient due to the two-photon continuum. The F coeffi cient 
is wavel e ngth- dependent and is taken from lAlleil (|l984l ) and 
iFerlandl (|l980l ). 



3 OVERVIEW OF SFR CALIBRATIONS AND 
PROPERTIES OF SFR TRACERS 

In order to compare our derived results with the previ- 
ous studies, in this part we will overview the previous re- 
sults about SFR calibrations in terms of luminosities of 
recombination- line, forbidden-line, UV and FIR continuum, 
and overview the properties (advantages and drawbacks) of 
these SFR tracers. 



3.1 Overview of SFR.vs.Luy 

The UV luminosity, Luv, is directly tied to the emission 
of the young SP. For convenience, it is often expressed as 
the linear relation with SFR by previous studies: SFR\jv — 
Fuv X Luv, where SFRuy and Luv are in units of Mq yr~^ 
and ergs~^ Hz~^. This kind of SFR calibration is valid only 
for galaxies with continuous SFR over timescales of 10* or 
longer (K98). Its advantage is that it can be applied to star- 
forming galaxies over a wide range of redshifts. Its drawback 
is that it is sensitive to extinction and the form of the IMF. 

Because different EPS models and methods are used, 
the obtained conversion factor Fuv between SFRw and 
Luv is different, being the difference as big as 0.3 dex (K98): 

• K98 has obtained Fuv = 1-4 x 10"^* by using the S55 
IM F with mass limits 0.1 and 100 Mq an d solar abundance. 

• iMadau. Pozzetti fc DickinsonI (|l998l . hereafter MPD98) 
have obtained Fuv = (1-25 x 10"2*,1.26 x 10"^**) at 
(1500 A, 2800 A) by using the S5 5 IMF , and Fuv = (2.86 x 
10~^^ 1.96 X 10"^*) by using the lScalol (|l986l ) IMF. In their 
studies the calibration factors are from models with an ex- 
pon entially decreas i ng SF R. 

• ICilbank et al.l (|2010l . hereafter GIO) have obtained 
Fuv = 0.71 X 10~^* by using the PEGASE code (assum- 
ing a constant SFR over 1 Gyr and a Kroupa IMF around 
solar metallicity). 



3.2 Overview of SFR.vs.Lua 

Also, the luminosity of Ha recombination line, Lua, is of- 
ten expressed as the linear relation with SFR: SFRua ~ 
Fhq X Lhq, where SFRna and Lua are in units of M0 yr^^ 
and erg s^^ . This equation traces the instantaneous SFR. Its 
primary advantage is its high sensitivity. In addition. Ha is 
typically so conspicuous that it can easily be detected. Its 
drawback is that it is sensitive to extinction, IMF and the 
assumption that all of the massive star formation is traced 
by the ionizing gas. The conversion factor Fua obtained by 
different studies is as follows: 

• K98 has obtained Fhc = 7.9 x 10""*^ by using the 
S55 IMF with a = 2.35 and the upper and lower mass 
limits Ml =0.1 and Mu = 100. M©. This value is ob- 
tained from a constant star formation model at 'equi- 
librium'. The similar value ha s also been obtained by 
iKennicutt. Tamblvn fc CongdonI (|l994l 'l and MPD98 when 
using the S55 I MF. This factor has been adopted by 
IShi. Gu fc P"en3 (|2006l ) and GIO. In the work of GIO, they 
have pointed out that the factor is similar to that obtained 
fror n the PEGASE code . 

• iBrinchmann et al] (|2004l , hereafter B04) have obtained 



6 F. Zhang, L. Li, Y. Zhang, X. Rang, Z. Han 



Fhq = 5.25 X lO"''^ by usi ns KOI IMF. This val ue is multi- 
plied by 1.5 in the works of lTremonti et aLll2004l and GIO in 
order to convert to the S55 IMF. A factor of 1.5 corresponds 
to the ratio of mass in the two IMFs for the same amount 
of i onizing radiatio n. 

• iMateus et~al] (j2007l ) have obtained Fhq = 5.22 x 
10~*^ X lO~° *'''~''f""'=', where r is the Petrosian magnitude 
representing the total galaxy flux, rgbre is the r— band fi- 
bre magnitude, and the last term is for aperture effect. This 
equation is derived by using the expression adapted from 
K98 for a Cha03 IMF (0 1 - lOOMp) and the prescriptions 
given bv lHoDkins etlil ( |200j. hereafter H 3). T his conver- 
sion factor has been used bv lHuang fc Gu| (|2009|). 

• K98 argued that when adopting the lScald (| 19861 ') IMF, 
the derived SFRua is approximately three times higher than 
that derived with a Salpeter IMF. 



3.3 Overview of SFJ?.vs.L[oii] 

The [OII]A3727A forbidden-line doublet is one of the 
strongest emission lines in the blue. It is accessible to optical 
observations over a wide range of redshifts. Therefore, it pro- 
vides a very us eful estimate of the SFR in distant galaxies 
(|Schaereilll999l . K98). Its disadvantage is less precise than 
that from Ha and may also be prone to systematic errors 
from extinction and variations in the diffuse gas fraction. 
Often, SFR is expressed as the linear relation with the lu- 
minosity L[oii]: SFR[ou] = F[oii] x I/[oii], where SFRioii] 
and L\ciu] are in units of M0 yr"'^ and ergs~^. 



^[oir 

K98 has obtained F[oii] 



(1.4 ±0.4) X lO""*^ by using 
L[ou]/LHa = 0.45. 

• H03 have obtained F[oii] ~ 3.37 x 10"''^ by using 
L[oii]/ Lua ~ 0.23 and the calibration factor Fua ~ 7.9 x 
10-*^ of K98. 

• GIO have obtained F[oii] ~ 3.95 x lO"''^ by assuming 
the ratio of extinguished [Oil] to Ha flux L[oii]/-f'Ha = 0.5, 
1 mag of extinction at Ha and the calibration factor Fh^ = 
7.9 x 10"*^ of K98. If neglecting the extinction, the GIO 



conversion factor F 



[Oil] 



1.57 X 10" 



3.4 Overview of SFR.vs.Lfir 

The FIR (10-300 /im) luminosity, Lfir, is only used to trace 
SFR of starburst galaxies with ages less than 10* years. Its 
drawback is that it i s also sensitive to IMF and e xtinction. 
Using the models of iLeitherer fc HeckmanI (| 19951 ) for con- 
tinuous bursts of age 10-100 Myr and the S55 IMF with the 
lower and upper mass limits Mi = 0.1 and Mu ~ 100 Mq, 
K98 has presented the conversion factor between SFRpm 
and Lfir: Ffir = 4.5 x 10"** Mq yr~ Verg 



rgs 



4 EFFECT OF BINARY INTERACTIONS ON 
SFR CALIBRATIONS AND COMPARISONS 

For the sake of clarity, we define eight sets of models. The 
two sets of Models A and B use the Yunnan EPS models with 
and without binary interactions, respectively. The six sets 
of Models G(Cr), D, E, F and G use the BC03, GISSEL98, 
PopSTAR, STARBURST99 and PEGASE models, respectively. 
In Table [2] we give the name of each set of models in the 1st 
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Figure 1. Relation between SFR and Lhq of E, SO, Sa, Sb, Sc 
and Sd galaxies (corresponding to r = 1,2,3,5,15 and 30 Gyr 
in Eq. |g] from top to bottom) for Models A-MS79 (black solid 
line-|-solid rectangles) and B-MS79 (red solid line-|-solid circles). 
The ages of galaxies are in the range of 1 Myr-15 Gyr. Also shown 
are the results of K98 (dashed line) and B04 (dotted line). 
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Figure 2. The evolution of Q(H) for Models A-MS79 (black 
linc+solid rectangles) and B-MS79 (red line-f-solid circles). From 
top to bottom are Burst, E, SO, Sa, Sb, Sc and Sd types of galax- 
ies. 



Table 2. Definitions of models. Bis and GR indicate binary in- 
teractions and gas recycle, respectively. 



Model 


EPS 


Bis 


IMF 


GR 


A 


Yunnan 


Y 


MS79/S55 


N 


B 


Yunnan 


N 


MS79/S55 


N 


C 


BC03 


N 


Cha03/S55 


N 


D 


GISSEL9S 


N 


MS79 


N 


E 


PopSTAR 


N 


K01/S55' 


N 


F 


STARBURST99 


N 


K937S55 


N 


G 


PEGASE 


N 


K937S55 


N 


Cr 


BC03 


N 


Cha03/S55 


Y 
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Figure 4. The distribution of stars in the logTcff-logg plane for the Burst galaxy type (SP) at an age of logt/yr =7.2. In eaeh panel 
black and red symbols arc for considering (i.e., Model A-X) and neglecting (i.e., Model B-X) binary interactions. Left and right panels 
are for Models A/B-MS79 and A/B-S55, respectively. 




Figure 6. Top panels arc the comparison in SFR-Lna relation between model (red symbols) and fitting (green symbols) results for 
Model A-MS79. Left and right panels are for using Eqs. I13l and ll4l respectively. Bottom panels are the residuals {logS F BAogS F Rai) as 
a function of logLna- 



column, the name of used EPS model in the 2nd column, the 
condition that binary interactions are taken into account in 
the 3rd column, the IMF in the 4th column, and the 5th 
column denotes the condition that the assumption of gas- 
recycle is used, where 'N' means e in Eq. [9] equals to zero 
(i.e., the gas could not be recycled into new star formation) 
and 'Y' means e — 1 (i.e., the gas could be recycled). In the 
seven sets of Models A-G e = 0, and in the set of Model Cr 
e = l. 

For each set of models two subsets are considered, de- 
pending on the IMF. To distinguish them, the name of used 
IMF (the 4th column) is the supplement to the name in the 
1st column of Table [51 

In this section we mainly discuss the effect of binary in- 
teractions on the SFR calibrations. Therefore, only Models 



A-MS79 and B-MS79 are used, which are based on the stan- 
dard Yunnan EPS models. The other models will be used in 
the next section to discussion the effects of IMF, gas-recycle 
assumption and EPS models on the results. 



In this section we obtain the number of ionizing photons 
(5(H), the luminosity of Ha recombination line Lua, the 
luminosity of [OIIJA3727 forbidden-line doublet L[oii], the 
UV fluxes at 1500 and 2800 A, Li,uv, and FIR flux Lpir of 
Burst, E, SO, Sa-Sd and Irr galaxies for both Models A-MS79 
and B-MS79. We also give the conversion factors between 
SFR and these SFR diagnostics. 
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Figure 7. Relations between SFR and uv of E, SO, Sa, Sb, Sc and Sd galaxies (corresponding to t = 1, 2, 3, 5, 15 and 30 Gyr in 
Eq. |9l from top to bottom) for Models A-MS79 (black solid line+solid rectangles) and B-MS79 (red solid line+solid circles). Left panel 
is for Li500i and right panel is for L2goo. Also shown are the results of K98 (grey dashed line), MPD98 (grey dotted line, open and solid 
triangles are for using the S55 and Scalo IMFs, respectively) and GIO (grey dash-dotted line). 
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Figure 8. Similar to Fig. [2] but the left panel is for L1500 and the right panel is for L 
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Figure 3. The stellar (solid line) and nebular (dotted line) spec- 
tra of the Burst galaxy type (SP) at an age of logt/yr =7.5 for 
Models A-MS79 (black+rectangles) and B-MS79 (red-|-circles) . 



4.1 SFR.vs.Lna 

In Fig. [1] we give the relation between log(5'-F'i?) and 
log(LHa) (note the logarithmic scale) of E, SO, Sa-Sc and 
Sd galaxy types in the range of 0.1 Myr-15 Gyr for Models 
A-MS79 and B-MS79. Also shown are the SFR(Lhc«) cali- 
brations of K98 and B04. If the aperture effec t is neglected, 
the conversion factor of iMateus et al.l |20o3) is similar to 
that of B04 (see Section [3]), so we do not show it in Fig. [1] 
for the sake of clarity. 

From Fig. [1] we see that the log(S'_Fi?) varies linearly 
with log(LHQ) within a certain age range (from logt/yr>7.5 
for E-type to logt/yr>8.5 for Sd-type in Model A-MS79, 
from logt/yr>7.1 for E-type to logt/yr>8.5 for Sd-type in 
Model B-MS79) for both Models A-MS79 and B-MS79. And, 
for an given log(LH<:«), the \og{SFR) of Model A-MS79 is 
smaller by an amount of '~-^0.2 dex than that of Model B- 
MS79 when \og{SFR) > - 11, i.e., the inclusion of binary 
interactions can lower the derived SFR in terms of Lhq. 
When comparing with the SFR of K98 and B04 at a given 
Lho, we find that the \og{SFR) increases from B04, K98, 
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Table 3. Conversion cocfBcicnts between SFR and Lhq for all 
models except for Models Cr-Cha03 and Cr-S55. 



LoglO(T^„) 

Figure 5. The evolution of a binary system (Afi = 13.97, M2 = 
5.44 M© and P = 28.0 day) from ZAMS to the end of heUum 
burning. Black and red symbols are for the primary and the sec- 
ondary, respectively. The red open circles denote the positions at 
an age of logt/yr = 7.2, corresponding to time in Fig. |4] 
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Figure 9. The Lpjji evolution of Irr galaxies (i.e., models with 
constant star formation, SFR=1M0) for Models A-MS79 (black 
solid line+solid rectangles), B-MS79 (red solid line+solid circles), 
A-S55 (green solid line+open rectangles) and B-S55 (blue solid 
line+open circles). Also shown is the result of K98 (green dashed 
line) . 



Model A-MS79 and Model B-MS79 in turn, and the discrep- 
ancy in the log{SFR) between K98 and Model A-MS79 is 
small. 

The reason that binary interactions lower the SFR at 
a given Lua is that they raise the number of ionizing pho- 
tons (3(H). In Fig. [5] we give the evolution of log(Q(H)) for 
Burst, E, SO-Sc and Sd galaxy types. From it we see that 
the log(Q(H)) of Model A-MS79 is significantly greater than 
that of Model B-MS79 for Bursts in the age range 6.7 < 
logt/yr < 8.4, the difference in log((5(H)) reaching ~2dex 
at an age of logt/yr=7.5. For E-SO-Sd galaxy types binary 
interactions also can raise log((5(H)) when logt/yr > 6.7. 

The raise of Q(H) is caused by the fact that binary in- 
teractions increment the UV spectrum at the corresponding 



Model 


Eq.[T3] 




Eg. [Til 






(Chq 


fHc) 






Ha' 


A-MS79 


-41.0556 


0.0096 


-41.4316 


1.0121 


0.0076 


B-MS79 


-40.8942 


0.0141 


-41.4904 


1.0193 


0.0104 


A-S55 


-40.7408 


0.0206 


-41.6298 


1.0289 


0.0151 


B-S55 


-40.6939 


0.0216 


-41.6182 


1.0301 


0.0159 


C-Cha03 


-41.3025 


0.0054 


-41.5413 


1.0077 


0.0040 


C-S55 


-41.0785 


0.0065 


-41.3640 


1.0093 


0.0048 


D-MS79 


-41.2936 


0.0055 


-41.5374 


1.0079 


0.0041 


E-KOl 


-40.7390 


0.0003 


-40.7296 


0.9997 


0.0002 


E-S55' 


-41.5047 


0.0002 


-41.4959 


0.9997 


0.0002 


F-K93' 


-40.4930 


0.0004 


-40.4773 


0.9995 


0.0004 


F-S55 


-40.8002 


0.0004 


-40.7850 


0.9995 


0.0003 


G-K93' 


-40.8519 


0.0653 


-38.7983 


0.9333 


0.0533 


G-S55 


-41.1342 


0.0724 


-38.8548 


0.9267 


0.0592 



Table 4. Conversion coefficients between SFR and L1500 (i.e., 
j = 1 in Eqs. 1161 and I17I I for all models except for Models Cr- 
Cha03, Cr-S55, E-KOl and E-S55'. 



Model 


Eq.[T6] 




Eq.[T7] 






(C1500 


CTISOO) 


(^'l500 


A' 




A-MS79 


-28.0681 


0.0219 


-28.4138 


1.0193 


0.0176 


B-MS79 


-28.0037 


0.0212 


-28.3104 


1.0171 


0.0178 


A-S55 


-27.8032 


0.0302 


-28.2842 


1.0272 


0.0242 


B-S55 


-27.7751 


0.0305 


-28.2403 


1.0263 


0.0249 


C-Cha03 


-28.1430 


0.0207 


-28.4825 


1.0193 


0.0154 


C-S55 


-27.9437 


0.0232 


-28.3225 


1.0218 


0.0172 


D-MS79 


-28.1410 


0.0208 


-28.4820 


1.0194 


0.0154 


F-K93' 


-27.7758 


0.0092 


-27.8092 


1.0019 


0.0091 


F-S55 


-27.8955 


0.0059 


-27.9120 


1.0009 


0.0059 


G-K93' 


-27.9137 


0.0635 


-27.7433 


0.9903 


0.0630 


G-S55 


-27.9968 


0.0413 


-27.6359 


0.9796 


0.0380 



ages. In Fig. [3] we give the stellar and nebular spectra of the 
Burst galaxy type (SP) at an age of logt/yr=7.5 for both 
Models A-MS79 and B-MS79. From it we see that binary 
interactions not only raise the stellar spectrum in the UV 
band, but also raise the nebular continuum (by an amount 
of ~2dex). 

The higher UV spectrum is caused by more hotter stars 
produced by binary interactions. In Fig. |4]we give the dis- 
tribution of stars in logrofr-log(? plane for Bursts (SPs) at 
an age of logt/yr = 7.2 for Models A-MS79 and B-MS79. 
From the left panel, we can see that binary interactions 
(i.e., Model A-MS79) indeed produce some hotter helium 
stars (logTcfT ~ 4.9 and log g ~ 5.8). These helium stars 
are originated from those binary systems, for which the 
primaries have evolved from the giant branch (GB, core- 
helium-burning) phase to the helium main-sequence (HeMS) 
phase due to the mass loss of the primaries. In Fig. [5] we 
give the evolution of such a binary system from zero age 
main sequence (ZAMS) to the end of helium burning in the 
logToff-logL plane. From it we see that the primary would 
evolve from GB to HeMS phase due to the mass loss, and 
the luminosity of the secondary increases rapidly due to the 
mass increase. In general, these binary systems have rela- 
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Table 5. Conversion coefficients between SFR and -L28OO (i-e., 
i = 2 in Eqs. 1161 and 1171 for all models except for Models Cr- 
Cha03, Cr-S55, E-KOl and E-S55'. 



Model 


Eq.[T6] 




Eq.[l7] 






(C28OO 


a-2800) 


(^2800 


A' 

^2800 


""2800) 


A-MS79 


-28.0227 


0.0674 


-29.2780 


1.0701 


0.0504 


B-MS79 


-27.9737 


0.0714 


-29.2873 


1.0735 


0.0540 


A-S55 


-27.7947 


0.0984 


-29.6921 


1.1073 


0.0729 


B-S55 


-27.7721 


0.1031 


-29.7416 


1.1115 


0.0773 


C-Cha03 


-28.0637 


0.0538 


-28.9987 


1.0534 


0.0390 


C-S55 


-27.8717 


0.0594 


-28.9053 


1.0597 


0.0428 


D-MS79 


-28.0625 


0.0539 


-28.9996 


1.0535 


0.0391 


F-K93' 


-27.7540 


0.0759 


-29.1760 


1.0805 


0.0569 


F-S55 


-27.8251 


0.0459 


-28.6379 


1.0458 


0.0351 


G-K93' 


-27.9478 


0.1704 


-30.5344 


1.1465 


0.1396 


G-S55 


-27.9483 


0.0963 


-29.1312 


1.0670 


0.0832 



lively large mass ratio q (>0.2) and small orbital period 

To quantitatively analyze the efTect of binary interac- 
tions on the calibration factor between SFR and Lna, we 
give a fitting relation between log(5'-Fi?) and log(Z/Hci ) when 
\og{SFR) > -11 and Alog{SFR) > 0.05 for Models A- 
MS79 and B-MS79 by the following form: 

SFRhci , Lua 



log 7 



= log- 



+ Cho 



(13) 



'(Moyr-i) ^(ergss-i) 

where SFRua means that it is calculated from Lua- The 
fitting coefficient (Chq) and the rms (cthq) are given in the 
2nd and 3rd columns of Table [3] respectively. This form of 
fitting (Eg. [T3ll is t he same as that of K98, MPD98 and 
iMateus et al. (|2007l ) . In fact, it is a linear relation between 
SFR and Lhc- 

From Fig. [T] we see that SFR does not vary linearly 
with Lua (it is easily seen when comparing the results of 
K98 and B04). Therefore, we also adopt the following form: 

SFRhc . I , I/Hq 



log 



(Mo yr-i) 



= Ana X log 



(ergs s 



+ Cho 



(14) 



The fitted coefficients (Ah^ & Chq) and the rms (cthc) are 
given in the 4th-6th columns of Table [3] In the work of H03, 
the relation between SFR and the U-band luminosity is fit- 
ted by the above form (see their Eq. 10). In the top panels of 
Fig.|6]we give the comparisons in the log(5'Fi?).vs.log(LHQ) 
relation between fitting (log SFR^t) and model data for 
Model A-MS79 in both cases of Eqs. [13] and [H In the 
bottom panels of Fig. [6] we give the residuals {logSFR- 
logSFRfit) as a function of logLua- 



4.2 SFR.vs.Li 



and SFR.vs.L2 



In Fig. [7] we give the relations between \og{SFR) and the 
logarithmic UV-luminosities at 1500 and 2800 A of E, SO-Sd 
types of galaxies for Models A-MS79 and B-MS79, also give 
the resuhs of K98, MPD98 and GIO. The resul ts of MPD98 
are obtained by using S55 with a = —2.35 and lScalcl (|l986l ) 
IMF. 

From the left and right panels of Fig. [71 we see that 
the differences in the SFR-Lii,oo and SFR-L2S0Q relations 
between Models A-MS79 and B-MS79 are small. By com- 
paring the SFR of K98, MPD98 and GIO at given L1500 



or 1/2800, we find that the \og{SFR) increase from GIO, 
Models A-MS79/B-MS79, MPD98-S55, K98 and MPD98- 
Scalo in turn (i.e., SFRqio > <S'F-Ra-ms79/b-ms79 > 
SFiiMPD98-S55 > SFRk_98 > •S'Fi?MPD98-Scaio), the dif- 
ferences between MPD98-S55 and K98 are small, and our 
resuhs lie between GIO and MPD98-S55. 

The reason that binary interactions do not affect SFR 
at given Z/1500 or L2800 can be seen from Fig. [8] in which 
we give the evolutions of I/1500 and Z/2800 of Burst, E, SO-Sc 
and Sd types of galaxies for Models A-MS79 and B-MS79. 
First, from the left panel of Fig. [8] we see that the L1500 of 
Model A-MS79 is greater than that of Model B-MS79 only 
for Bursts at ages of logt/yr ~ 9 (the maximal difference 
is ~1 dex at an age of IGyr). For the other galaxy types 
the difference in the Lisoo is small. The reason for this is 
that the contribution of the SP with an age of ~lGyr to 
the ISEDs is smaller for E, SO-Sd galaxies according to 

F{t)^ [ i,{t-t')fsp{t')At', (15) 
Jo 

where F{t) is the galaxy spectrum at time t, ^{t — t') is 
the SFR at t -t' (see Eq. ^ and /sp(t') is the fiux of SP 
with an age of t' . In addition, from the right panel we see 
that the difference in the L2800 between Models A-MS79 
and B-MS79 is insignificant for all galaxy types at all ages. 
Therefore, binary interactions would not affect significantly 
the SFR-L1500 and SFR-L2S00 relations. 

Also we fit the log{SFR)-log{Li50o) and log{SFR)- 
log(Z/280o) relations by using the expressions 



log 



SFR,, 



(Ms yr-i) 



= log 



(ergs s~^ Hz" 



and 



log 



SFR,i 



(Mo yr- 



Ai,uv X log 



(ergs s~i Hz~^ 



(16) 



+ C:,uv, (17) 



where i=l denotes the wavelength A = 1500 A, i=2 means 
A — 2800 A, and SFRi^jjy means that it is from the i-th UV 
luminosity Li.uv. The fitting coefficients between log{SFR) 
and log(Li5oo) (C,,uv, C'i_uv, A'^ uv) and rms {ui,vv,r^'i,vv) 
are given in Table O and those between \og{SFR) and 
log(L280o) are given in Table [5] 

4.3 SFR.vs.L[oni 

The luminosity of the [OII]A3727A forbidden-line doublet, 
L[oii], is indirectly related to the ionizing luminosity. Often 
it is obtained by using the ratio of L[oii]/-f'Ha, which is ob- 
tained empirically. In the work of K98 L[oii]/iHci = 0.45, 
in the work of H03 it is 0.23, in the work of GIO the value 
is 0.5, and in the pegase code the value is 3.01/2.915. In 
this work we use the value of 0.23 used by H03 to obtain 
L[oii]. Because the fixed L[ojj^/Luct ratio is used, the plot 
of the 5'-F'ii.vs.L[oii] is similar to Fig. [1] At a given SFR, 
the log(L[oii]) is smaller than log(LHa) by an amount of 
log(l/0.23), i.e., the calibration curve moves upwards by 
log(l/0.23). Therefore we do not give the plot of SFR-L[ou] 
relation. Combining the conclusion made in Section [4.11 we 
know that binary interactions make the log{SFR) in terms 
of I/[oii] larger by about 0.2 dex at a given I/[oii]- 

Also, we give a fitting relation between SFR and L[oii] 
by the following two forms: 
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Figure 10. Similar to Fig. [T] but including the results of Models 
A-S55 (green solid line+open rectangles) and B-S55 (blue solid 
line+open circles). 
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Because the fixed L^oii]/Lua ratio is used, the fitting 
coeflicient Ajon] = Aljj^, C[oii] = Chq - log(0.23) and 
C(oii] = Ch^ - log(0.23) (cf. Eqs. \lMM- The actual num- 
bers can immediately be computed from the values in Table 
3. 



4.4 SFR.vs.Lfir 

In the computations of the above three SFR diagnostics 
(Section l4.H4.3[) . the calibration factors are from those mod- 
els with an exponentially decreasing SFR (i.e., Eq. |9]l, while 
the calibration between SFR and FIR luminosity is from 
the models with constant SFR under the assumption of the 
bolometric luminosity Lbol ~ ipia- 

In Fig. [9] we give the Lpir evolution of Irr galaxies 
(i.e., models with constant SFR) for Models A-1VIS79 and 
B-MS79. Also shown are the results of K98. In all cases the 
total mass of the model galaxy is normalized to 1 Mq . From 
this figure we see that the difference in the Lpir evolution 
between Models A-MS79 and B-MS79 is small, that is to 
say, binary interactions almost do not vary the conversion 
factor between SFR and Lfir- 



5 OTHER FACTORS ON SFR CALIBRATIONS 

The calibrations between SFR and Lua, iisoo, ^2800, i[oii] 
and LpiR, which have been given in Section |4l can be af- 
fected by the adoption of different EPS models, IMF and 
the assumption of gas recycle. In this section, we will dis- 
cuss the effects of these factors on the above mentioned SFR 
calibrations by introducing the Models A-S55, B-S55, G- 
K93' and G-S55. For these models, the Lua, iisoo, -f'2800, 
L[oii] and Lfir are calculated, and the calibration factors 
are presented in Tables. |3l |4] and [5l respectively. 
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Figure 11. Comparison of the IMF between IVIodels A/B-MS79 
(rod, solid lino) and A/B-S55 (cyan, dashed line). The thick lines 
arc for the mass of binary system, the upper and lower lines are 
for the masses of the primary and secondary, respectively. Note 
in this figure logAf is obtained by using the EFT's approximation 
to the MS79 and S55 IMFs. 
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Figure 12. Similar to Fig. |2] but also shown are the results 
of Models A-S55 (green line-|-open rectangles) and B-S55 (blue 
line-(-opcn circles). 

5.1 The effect of IMF 

To analyse the effect of IMF on these calibrations, the results 
of Models A-S55 and B-S55 are needed to be combined with 
those of Models A-MS79 and B-MS79 (Section H]). Models 
A-S55 and B-S55 are built by using the Yunnan models with 
the S55 IMF and the assumption that the masses of the two 
component stars in a binary system are uncorrelated. 

5.1.1 IMF on SFR. vs. Lhc 

In Fig. [To] we give the comparison in the SFR-Lua relation 
among Models A-MS79, A-S55, B-MS79 and B-S55. Also 
shown are the results of K98 and B04. 

By comparing the SFR-Lna relation between Models 
A-MS79 and A-S55, and that between Models B-MS79 and 
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Figure 14. Similar to Fig. [T] but including the results of Models A-S55 (green solid line+open rectangles) and B-S55 (blue solid 
line+open circles). 
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Figure 15. Similar to Fig. [8] but including the results of Models A-S55 (green line+open rectangles) and B-S55 (blue line+open circles). 



B-S55, we see that the log(I/Ha) of Model A-S55 is smaller 
than that of Model A-MS79 by an amount of 0.4 dex, and 
that the prediction for Model B-S55 is smaller by ~0.2 dex 
than that of Model B-MS79 for a given log{SFR). This is 
partly caused by the fact that Models A/B-S55 produce less 
massive-stars (see Fig. 1111) and less Q(H) (see Fig. I12|) . In 
Fig. [11] we present a comparison of the IMF between Models 
A/B-MS79 and A/B-S55, where the S55 and MS79 IMFs 
are obtained by using the EFT's approximation (i.e., Eqs.[T] 
and |4]) and are different from those given by Eqs. [2] and [5] 
In Fig.[T2]we present the Q(H) evolution of all galaxy types 
for Models A-MS79, B-MS79, A-S55 and B-S55, respectively. 
From them we indeed see that the number of massive stars 
and therefore (3(H) at early ages of Bursts for Models A/B- 
S55 is less than those for Models A/B-MS79. 

Moreover, by comparing the result between Models A- 
S55 and B-S55 we find that the difference in the SFR-Lna 
relation between them is small, which is significantly dif- 
ferent from that between Models A-MS79 and B-MS79. The 
difference in the SFR-Lhc relation between Models A-MS79 
and B-MS79 is caused by the difference in the (5(H) of Bursts 
in the age range 6.7 < logt/yr < 8.4, while the difference 



in the (5(H) between Models A-S55 and B-S55 is small for 
Bursts at all ages (see Fig. [T2I . Why Model A-S55 could not 
produce more Q(H) than Model B-S55, like Model A-MS79? 

This is because less hotter helium stars (logTefi 4.9 
and logg ~ 5.8) could be produced in the range 6.7 < 
logt/yr < 8.4 for Model A-S55 (see the right panel of Fig.|f| 
and the comparison with the corresponding left-right panel) . 
We can understand it from the discussion in Section |3?T] and 
from Figs. Illland ll3l First, from the discussion in the Sec- 
tion 14.11 we know that those hotter helium stars (present 
in the left panel of Fig. |3| evolve from those initial bi- 
nary systems with the relatively large primary-mass (Mi), 
large mass-ratio (q) and small orbital separation. Then, from 
Fig. [TT]we know that the number of massive primary-stars 
becomes to be less for Model A-S55. Furthermore, from 
Fig. [13] we know that the number of binary systems with 
relatively large AIi and large q becomes to be less for Model 
A-S55. In Fig. [I3]we give the the number of binary systems 
in the primary-mass range A/i ^ Mi + dMi and the mass- 
ratio range q ^ q + dqfor Models A/B-MS79 and A/B-S55. 
The results of using the S55 IMF and assuming that the 
masses of the two component stars are correlated (the same 
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Figure 13. The number of binary systems in the ranges Mi — > 
A'/i + dA4i and q q + dq (dg=0.02) per million binary systems 
(dAf) for Models A/B-MS79 (red) and A/B-S55 (cyan). The grey 
circles denote the grid (Mi, q), and the distance from grid to the 
right point denotes the number of binary systems. For compari- 
son, we also give the results of using the S55 IMF and the assump- 
tion that the masses of the two component stars are correlated 
(blue). 

as that in Models A/B-MS79) are also represented. For the 
sake of clarity, in Fig. [13] we only give the number of binaries 
with the primary-mass 8 < Mi < 45M0 . These binaries are 
likely to evolve into hotter helium stars in the range 6.7 < 
logf/yr < 8.4, and the MS lifetimes of single stars with the 
upper and lower masses (8 & 45 Mq) approximately corre- 
spond to the ages of logf/yr=6.7 and 8.4 (the MS lifetime 
of M = 3.5Mo is logt=8.4, the MS lifetime of star with 
M = 38. M0 is \ogt—6.7 at solar metallicity). From it we see 
that the number of binary systems with larger Mi and q 
is less for Model A-S55. Thus Model A-S55 could not pro- 
duce more hotter helium stars and Q(H). Moreover, from 
Fig. [13] we also see that the number of binaries with larger 
Ml and q of Model A-MS79 is similar to that by using S55 
IMF and assuming that the masses of the two component 
stars are correlated. Therefore, the real reason why Model A- 
S55 could not produce hotter helium stars is the assumption 
about the masses of two component stars in binary systems. 

From the above discussion, we conclude that the 
SFR.vs.Lua calibration is affected not only by the adop- 
tion of the different IMF, but also by the assumption about 
the masses of the component stars in binary systems. 



5.1.2 IMF on SFR. vs. Li5oo and L2S00 

In Fig. [TJwe give the comparisons in the SFR-L1500 and 
SFR-L2S00 relations among Models A-MS79, A-S55, B- 
MS79 and B-S55. Also shown are the results of K98, MPD98 
and GIO. 

By comparisons, we find that the log(Li5oo) and 
log(Z/280o) of Models A/B-S55 are smaller than the corre- 
sponding ones of Models A/B-MS79 by about 0.2 dex at a 
given \og{SFR), the conversion curves move upwards from 
the lines of K98 and MPD98-S55. For the reason we pro- 
vided in Section 15.1.11 (less massive stars), and from the 



results displayed in Fig. [15] (which represents the UV lu- 
minosities for all galaxy types predicted by the Models A- 
MS79, B-MS79, A-S55 and B-S55). We see that the L1500 
and 1/2800 of Models A/B-S55 are smaller than the corre- 
sponding ones for Models A/B-MS79 at all ages for Bursts, 
leading to smaller I/1500 and I/2800 for E, SO-Sd galaxies and 
therefore larger conversion factors between SFR and Lisoo 
and between SFR and 1/2800- 

From the left panel of Fig. 1141 we see that the difference 
in the SFR-L1500 relation between Models A-S55 and B-S55 
is small, which is similar to the case of using MS79 IMF 
(see Fig. [7]). The reason is the same given in the discussion 
presented in the 3rd paragraph of Section r4.2l i.e., the I/1500 
of Models A-S55 is smaller than that of B-S55 only for Bursts 
in the age range 8.75 < logf/yr < 9.2, and the maximal 
difference between them is ~1 dex (also see the left panel of 
Fig. I15p . This would not affect significantly the L1500 of E, 
SO-Sd galaxies. Therefore IMF would not affect significantly 
the SFR.vs.Lisoo conversion. 

Moreover, from the right panel of Fig. 1141 we see that 
the difference in the SFR-L2S00 relation between Models A- 
S55 and B-S55 is also small. This is because that IMF does 
not affect the Z/2800 for Bursts at all ages. 



5.1.3 IMF on SFR.vs.Lioii] 

Because we use the fixed L^oiii/Lua ratio, the effect of IMF 
on the SFR.vs.L^oii] conversion is the same as that on the 
SFR.vs.Lua conversion, i.e., the conversion factor between 
SFR and Z/[oii] of Model A-S55 is larger than that of Model 
A-MS79 by ~0.4dex, the conversion factor of Model B-S55 
is larger than that of Model B-MS79 by ~0.2 dex, and the 
difference in the SFiJ.vs.Ljoii] conversion between Models 
A-S55 and B-S55 is small. 



5.1.4 IMF on SFR.vs.Lfik 

In Fig. [9] we give the Lfir evolution of Irr galaxies for Mod- 
els A-S55 and B-S55. From it we see that the Lfir of Mod- 
els A/B-S55 is lower than that of Models A/B-MS79 by 
an amount of ~0.2 dex. Therefore, the conversion factor be- 
tween SFR and Lfir of Models A/B-S55 is larger than that 
of Models A/B-MS79 by an amount of ~0.2dex. 



5.2 The effect of gas-recycle assumption 

In this part we use the Models C and Cr to discuss the 
effect of gas-recycle assumption on these SFR calibrations. 
In Models C and Cr the value of e (see Eq. O is set to and 
1, respectively, i.e., in Model C the gas can not be recycled 
into new star formation, while in Model Cr it can. 

In Fig. [TS] we give the calibration between 5*^7? and 
Lhc. for Models C-Cha03, C-S55, Cr-Cha03 and Cr-S55 (left 
and right are for Models C and Cr, respectively). By com- 
parison, we see that \og{SFR) does not vary linearly with 
log(LHa) within a certain SFR (age) range for E, SO-Sd types 
of galaxies when considering the gas-recycle assumption (i.e. 
Models Cr-Cha03/S55), and that log(Z/Hc«) of Models Cr- 
Cha03/S55 is greater than that of Models C-Cha03/S55 at 
a given log(5'Fi?) (easily seen by comparing with the line 
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Figure 17. Similar to Fig. [2] but for Models C-Cha03 (black 
line+solid rectangles) and Cr-Ch03 (red line+solid circles) and 
only for Sb-type galaxies. 



of B04). This is because the inclusion of the gas-recycle as- 
sumption produces more (5(H) and larger Lua- In Fig. [17] 
we give the evolution of (5(H) of Sd-type galaxies for Models 
C-Cha03 and Cr-Cha03. 

Moreover, from Fig. [16] we see that the difference in 
the log(S'-F-R)-log(LHQ) relation between Models C and Cr 
increases with decreasing SFR (increasing age). For exam- 
ple, the line of Model C-Cha03 overlaps that of B04 at 
log(5F_R) ~ -11, while that of Cr-Cha03 is shifted to the 
right of the line of B04 (~0.1 dex). The reason for this is that 
the difference in the (5(H) increases with age (see Fig. I17|) . 

For the Lisoo and I/2800, the situation is similar to that 
of Z/Ha- For the sake of size, we do not give the plots for 
them in this paper. The inclusion of gas-recycle assumption 
can lower the derived SFR in terms of Li,uv by an amount 
of 0.05 dex at log{SFR) ~ -12. 

For the -L[oii], the situation is the same as that of Lua- 
As for the Lfir, the inclusion of gas-recycle assumption al- 
most does not affect the conversion factor between SFR and 
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Figure 19. Comparison of IMF (including S55, S55', K93', KOI 
and MS79). The detailed descriptions are presented in Section 2. 



5.3 The effect of EPS models 

We use all Models, except for the set of Model Cr, to analyse 
the adoption of different EPS models on these SFR calibra- 
tions. These models use the EPS models of Yunnan, GISSEL98, 
BC03, PopSTAR, STARBURST99 and PEGASE, respectively. 

To decrease the influence of IMF, we divide the com- 
parisons into two groups. One group uses the S55 IMF and 
another group uses NON-S55 (including KOI, K93', MS79 
and Cha03) IMF. Moreover, we only compare our results 
with those studies by using the same (similar) IMF. 

5.3.1 EPS models on SFR.vs.Luc 

In Fig. [18] we give the comparison in the SFR-Lua relation 
for all models by using S55 and NON-S55 IMFs. 

At first, from the left panel of Fig.[TS]we see that Models 
A-S55, B-S55 (Yunnan) and F-S55 {starburstqq) give the 
similar SFR-Lua calibration and the corresponding curves 
lie above the line of K98. The calibration curve of Model 
E-S55' (PopSTAR) locates below the line of K98 and is the 
lowest one; the calibration curve of Model C-S55 (bgos) is 
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Figure 18. Similar to Fig. [T] Left panel is for all Models using S55 IMF (including A-S55, B-S55, C-S55, E-S55'; 
Right panel is for all Models using NON-S55 IMFs (including A-MS79, B-MS79, C-Cha03, D-MS79, E-KOl, F-K93^ 
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F-S55 and G-S55). 
and G-K93'). 
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Figure 20. Similar to Fig. [7] Left and right panels correspond to L1500 and 1/2800 1 respectively. Top panels represent all the models 
using S55 IMF (including A-S55, B-S55, C-S55, E-S55', F-S55 and G-S55). Bottom panels show all the models using NON-S55 IMFs 
(including A-MS79, B-MS79, C-Cha03, D-MS79, E-KOl, F-K93' and G-K93'). 



close to that of K98; and the calibration curve of Model 
G-S55 (pegase) lie between Models C-S55 and E-S55' and 
does not display a linear calibration relation. The difference 
in the SFR.YS.Lua calibration caused by the adoption of 



different EPS models can reach ~ 0.7 dex when using the 
S55 IMF. 

From the right panel of Fig. [TH] we see that the cali- 
bration curves of Models C-Cha03 (bcos) together with D- 
MS79 {GISSEL98) overlap the line of B04, and the calibration 
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Figure 21. Similar to Fig. |9] Left panel represents all the models using S55 IMF (including A-S55, B-S55, C-S55, E-S55', F-S55 and 
G-S55). Right panel shows all the models using NON-S55 IMFs (including A-MS79, B-MS79, C-Cha03, D-MS79, E-KOl, F-K93' and 
G-K93'). 



curves of the other models locate above that of B04 (from 
bottom to up, A-MS79, B-MS79, E-KOl. F-K93'). The dif- 
ference in the SFR.vs.Lua calibration is ~ 0.9 dex when 
using NON-S55 IMF. This difference is partly caused by the 
difference in the IMF. In Fig. [19] we give the comparison 
among the S55, S55', KOf , K93' and MS79 IMFs. 



5.3.2 EPS models on SFR.vs.Lisoo and SFR.VS.L2S00 

In Fig. [20] we give the comparisons in the SFR-L1500 and 
SFR-L2S00 relations for all models except for the set of 
Model E (popSTAR), for which the ISEDs are not provided. 
Top and bottom panels are for the results based on S55 and 
NON-S55 (including KOI, K93', MS79 and ChaOS) IMFs, 
respectively. 

At first, from the top panels of Fig. 1201 we see that 
the SFR.VS.L1500 and SFR.VS.L2S00 conversion factors of 
all models with the S55 IMF (A-S55, B-S55, C-S55, F- 
S55 and G-S55) are similar to those of K98 and MPD98- 
S55 (i.e., the resuhs with the S55 IMF). The differences 
in the SFR.vs.Li^oo and SFR.VS.L2800 conversion factors, 
caused by the adoption of different EPS models, are less 
than ~0.3 dex when using the S55 IMF. 

From the bottom panels of Fig. 1201 we see that all 
models with NON-S55 IMF (A-MS79, B-MS79, C-Cha03, 
D-MS79, F-K93' and G-K93') give small SFR-Li^oo and 
SFR-L2S00 conversion factors than the corresponding ones 
of MPD98-Scalo and similar values to those of GIO (for 
MPD98-Scalo and GIO, the NON-S55 IMF is used). 

The differences in the SFR.vs.Li^oo and SFR.YS.L2S00 
conversion factors, caused by the adoption of different EPS 
models, can reach ~0.2 dex when using NON-S55 IMF. This 
is partly caused by the difference in the IMF. 



5.3.3 EPS models on SFR.vs.Lioii] 

The effect of EPS models on the SFR-L[oii] relation is the 
same as that on the SFR-Lhc relation, i.e., the difference 
in the 5'-F'-R.vs.L[oii] conversion factor can reach ~ 0.7 dex 



when using the S55 IMF, and ~ 0.9 dex when using the 
NON-S55 IMF. 



5.3.4 EPS models on SFR.vs.Lfir 

In Fig. [21] we give the Lpir evolution of Irr-type galax- 
ies (i.e., models with const SFR) when using the S55 and 
NON-S55 IMFs for aU models except for the set of Model 
E (PopSTAR). From this figure we see that the difference in 
the EPS models can cause the difference of 0.4 dex in the 
SFR.vs.Lfik conversion factor when using the S55 IMF, 
and the difference of 0.8 dex when using the NON-S55 IMF. 



6 SUMMARY 

We use the Yunnan EPS models with and without binary 
interactions to present the luminosities of the Ha recombi- 
nation line, the [OII]A3727 forbidden-line doublet, the UV 
(at 1500 and 2800 A) and the FIR continuum for Burst, E, 
SO, Sa-Sd and Irr galaxies, and present the calibrations of 
SFR in terms of these diagnostics. 

By comparison, we find that binary interactions lower 
the SFR.vs.Lhq and SFR.vs.L[oii] conversion factors by 
~0.2dex, and do not significantly vary the SFR.vs.Li,uv 
(at 1500 and 2800 A) and SFR.vs.Lpir calibrations. 

We also consider the effects of IMF, the gas-recycle as- 
sumption and EPS models on these calibrations. By compar- 
ison, we find that the SFR.vs.I/Ha and SFR.vs.L[oii] con- 
version factors of Models A/B-S55 are larger by 0.4 and 
0.2 dex than the corresponding ones of Models A/B-MS79, 
and that the SFR.vs.Z/i,uv and SFR.vs.Lpir conversion fac- 
tors are larger by 0.2 dex. By comparing the results between 
Models C and Cr, we find that the inclusion of gas-recycle 
assumption only lowers the SFR calibrations at faint SFR. 
Also we use the other EPS models {bC03, GISSEL9S, Pop- 
STAR, PEGASE and STARBURST99) to obtain these SFR cal- 
ibrations. By comparison, we find that the differences in 
the SFR(LHa) and SFR(L[oii]) calibrations reach ~ 0.7 and 
0.9 dex, the difference in the SFR(_Lfir) calibration reaches 
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0.4 and 0.8 dex, and the differences in the SFR(I/i,uv) cali- 
bration reach 0.3 and 0.2 dex when using S55 and NON-S55 
(partly caused by the difference in the IMF) IMFs, respec- 
tively. 

At last, in this paper we give the conversion coefficients 
between SFR and these diagnostics for all models. In this 
paper we have only considered the effects of binary inter- 
actions for solar metallicity galaxies - more detailed studies 
will be given. 
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